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• • Abstract 

^ ' Let w{t,x) := {u,v){t,x), t > 0, x £ M", be the M^-valued spatial-temporal 

random field w = (n, v) arising from a certain two-equation system of fractional 
, kinetic equations of reaction-diffusion type, with given random initial data n(0, x) 

' and v{0,x). The space-fractional derivative is characterized by the composition 

^] I of the inverses of the Riesz potential and the Bessel potential. We discuss two 

Tj" ■ scaling limits, the macro and the micro, for the homogenization of w(t,x), and 

o ; 

, prove that the rescaled limit is a singular field of multiple Ito- Wiener integral 

type, subject to suitable assumptions on the random initial conditions. In the 
^ . two scaling procedures, the Riesz and the Bessel parameters play distinctive roles. 

' Moreover, since the component fields u, v are dependent on the interactions present 

within the system, we employ a certain stochastic decoupling method to tackle this 
components dependence. The time-fractional system is also considered, in which 
the Mittag-Leffler function is used. 
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1 Introduction 

The purpose of this paper is to present a certain homogenization theory associated with 
the following linear fractional kinetic system of reaction-diffusion type 



in the above, /ij > 0, < /3 < 1, < a < 2, < 7, and t > 0, x G W. The parameter j3 
denotes the time-fractional index, and a, 7 denote the space-fractional indices, for which 
we refer as the Riesz parameter and the Bessel parameter respectively (see [281 V.l and 




When a, 7) = (1, 2, 0), the system (11. ip is reduced to a classical reaction-diffusion 
system. The time-fractional index /3 < 1 means sub-diffusive (super-diffusive in case 
/3 > 1, which we do not study in this paper; see the remark in Section 6). The spatial- 
fractional Riesz index a means the jumps of the evolution, and Bessel index 7 means 
the tempering of large jumps; see the now-classic book of Stein |28l V.l and V.3] for 
precise mathematical explanations. To our knowledge, fractional kinetic equations of 
Riesz-Bessel type appear firstly in Anh and Leonenko |21I2]; abundant subsequent works 
in this direction by the authors and collaborators can be seen in [H El El [3 [TH [17] 
and the references therein. The two papers [71 [T7] with external potentials are partic- 
ularly related to the study of this paper. We should mention that fractional operators 
with two fractional parameters are natural mathematical objects to describe long-range 
dependence and/or intermittency; one can find data exhibiting such characteristics in 




(1.1) 




V.3]). 
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a large number of fields including economics, finance, telecommunications, turbulence, 
and hydrology. 

In this paper, we consider the system (11. ip with /ii = /i2 = /i > and with the 
random initial data Uq and Vq, of which are independent and each one has a certain 
long-range dependence in its random structure. This paper is along the jH] on single 
fractional kinetic equation; yet our results in this paper are with the following novel 
features. Firstly and most importantly, we study two scaling procedures, the macro and 
the micro, of the homogenization of the associated spatial-temporal random solution- 
field, in which the Riesz and the Bessel parameters play distinctive roles; our result 
on the micro-scaling is new, even for the single equation case; this micro-scaling makes 
use of both the Riesz and the Bessel indices and also needs the rescaling on the initial 
data. We feel that our result in this micro-scaling may capture the proclaimed feature 
of the intermittency of the random motions. Moreover, due to the interactions present 
within the system, the components fields u, v are dependent (even we have assumed the 
independence of the initial data Uq^Vq), and we employ a certain stochastic decoupling 
method to tackle this components dependence. Our study may show how the theory of 
Riesz potentials and Bessel potentials, as in the Chapter V of Stein [28], may appear 
significantly in the homogenization of random fields. 

The study on the single P.D.E. with random initial condition can be traced back to 
[T5] and [26], and then has a long active development; we refer to the citations in the 
above and the references therein. There also has very significant progress on Burgers' 
equation with random initial data; see the monograph of Woyczyhski [3U] and the Chap- 
ter 6 of Bertoin [H]. Whilst, to our knowledge, relevant study on P.D.E. system with 
random initial data seems few in previous literatures, except the works of Leonenko and 
Woyczyhski [191 [20] multi-dimensional Burgers' random fields (Burgers' Turbulence). 

The results of this paper show that, for the fractional kinetic system (II. ip with 
suitable random initial data, the rescaled field w^(t, x) in the micro-scaling (which we 
mean et and e J, 0) of homogenization both the Riesz parameter a and the Bessel 
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parameter 7 play their roles; while in the macro- scaling (which we mean ^ and e I 0) 
only the Riesz parameter a plays the role. Nevertheless, in either case the limiting field 
is a singular field of multiple Ito- Wiener integral type. Furthermore, the component 
fields u, V are dependent, due to the interactions present within the system, we employ a 
certain stochastic decoupling method to tackle this components dependence. The method 
itself could be potentially important in the future study on some random systems, for 
example the gradient system of Hamilton- Jacobi equation with random initial data, as 
in pm p. 173] ; we notice that the decoupling has been traditionally used in solving 
differential equation systems. 

The underlying idea in this paper is motivated by those works in [21 |3l [71 [T71 [18] 
and the references therein. Namely, we use the spectral representations to describe the 
sample field arising from the initial data, and the relations between Hermite polynomials 
and homogeneous chaos associated with the initial data, to get representations for the 
limit field in terms of multiple Ito- Wiener integrals. From limit theorems point-of-view, 
our results, also those in the above citations, belong to the realm of non-central limit 
theorems for convolution type integrals, in which the papers [221 [12] are pioneering; see 
also the monograph of Major[25] and survey papers in the special volume edited by 
Doukhan, Oppenheim and Taqqu [10] 

The paper is organized as follows. In Section 2 we give the explicit solution of the 
system (11.11) . In Section 3 the initial data are assumed to be stationary random fields, 
and we discuss the covariance structure of the resulting solution-vector random field of 
(11.11) . subject to the specified random initial condition; we show that the spectral method 
is suitable in describing our random field relative to the space-time parameter. We also 
introduce the initial field to be a certain subordinated Gaussian random field generated 
by a class of non-random functions whose variables are relative to the spatial parameter. 
In the main Sections 4 and 5, we consider (II. ip with the usual time-derivative and the 
fractional spatial-derivative characterized by the Riesz and the Bessel parameters. We 
present the homogenization of micro-scaling in Section 4 and the less subtle macro- 



4 



scaling in section 5, respectively. In Section 6 we provide extensions of the results in 
Sections 4 and 5 to the time- fractional /3 < 1, in which we need to use the Mittag-Leffer 
function. The proofs of all our results are given in Section 7. 

Acknowledgement. The authors are grateful to the inspiring lectures of Professor 
W.A. Woyczyhski at National Taiwan University for the perspective on Mathematical 
Theory of Fractional P.D.E. 



2 Preliminaries 

To begin with, we rewrite the system ( II. ip . with /ii = /i2 = /i > 0, a > 0, 7 > and 
/3 = 1, in the matrix form as follows: 

subject to some initial conditions 



n(0,x) \ _ f Uo{x) 



X e M", (2.2) 



where u = u(t,x), v = v{t,x), t > 0, x E M", A is the n- dimensional Laplacian, and B 
is a 2 X 2 matrix. 

The Green function G{t, x; a, 7) associated with the operator dt + /u(/ — A) 2 (— A)t 
is represented via the spatial Fourier transform as follows; see, f2E\ Chapter 5] or [51 
Section 2]. 

[ e'<^'^>G(t,x;a,7)(ix = exp[-/it|A|"(l + lA^i], A e M", (2.3) 

where < ■, ■ > denotes the inner product on R". 

In order to get a explicit representation for the solution of fl2.ip . we impose the 
following assumption on the matrix B. 

Condition A. Suppose the matrix [%]i<i.j<2 is diagonalizable, i.e., the matrix B can 
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be written as 

B=(^'' ] = PDP-\ with P={^'' \ , (2.4) 
\ hi b22 J \ P21 P22 J 

where P is a real-valued non-degenerate eigenvector matrix associated with the matrix 
B, and D = diag{di,d2), di, ^2 G M, where dj is the eigenvalue associated with the 
eigenvector {pi.j,p2.j)'^ (here and henceforth, T denotes the transpose). Without loss of 
generality, we suppose that det{P) = 1. 

Under Condition A, the Cauchy problem (12. II) (12.21) has the unique solution given by 

^ - Q{t; d,, d^) { 1 , t > 0, X G R^ (2.5) 



v{t^x) I \ V{t,x) 



where 



Qit;d,,d2):=P\ ^ \P^\ (2.6) 



e' 



and U{t,x), V{t,x) are determined by 



Mi^ X = / G{t,y;a,i) . \dy, (2.7) 

y\t,x) j jRn y y^i^x-y) j 

where the Green function Git^y] 0,7) is defined in (12. 3p . 

For completeness, we give the sketchy proofs of (12.51) . Firstly, by taking the spatial 

Fourier transform on both sides of (12.11) . under Condition A, we have 

^ (i, A) = (-/i|Ar(l + |A|2)U B) I ! 1 (t. A) 



dt 



V 



' -,i\\\'{\ + \\\^)i+<h] \v ' ^ ' 



Thus, 

(t,A) 



u 

V 



, ^ , -w|A|°(l + lARi \ ^ ^^ u , , 

) ( 1(0. A). (2. 
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Finally, f l2.5p and f l2.7p are followed by taking the inverse Fourier transform on both sides 
of fl2.8p and using the representation fl2.3p . Additionally, by fl2.3p we can also observe 
that 

G{t, x; a, 7) dx = 1, for any t > 0. (2.9) 

3 Correlated random structures 
3.1 general random initial data 

Firstly, we set {Q, J-", V) to be an underlying probability space, such that all random 
element appeared in this paper are measurable with respect to it. 
The following condition is imposed on the initials, in which and henceforth. 
Condition B. Let Uo{x) = ?7i(x) = rii{x,u) and Vo{x) = 772(2:) = ri2{x,u), x G 
1^*^,0; G Q, be independent, and distributed as two real, mean-square continuous, ho- 
mogeneous and isotropic random fields defined on the underlying complete probability 
space {fl,J-',V). We assume that Erij{x) = 0, Var(?7j(x)) = 1, and covariance functions 

i?,^.(x)=^,^(|x|):=Cov(r/,(0),r/,(x))= / e'<^'^> F,{dX), jG{1,2}, 

where the last equality is guaranteed by Bochner-Khintchine theorem and Fj{-) is the 
spectral measure corresponding to the field ?7j(-) for j G {1,2}, respectively. 

Under Condition B, in view of Karhunen's Theorem (see, for example, Gihman and 
Skorokhod [13], pp. 208-230), there exist complex- valued orthogonally scattered random 
measures Zpj, j G {1,2}, such that the random fields rjjix), j G {1,2}, have the 
following spectral representations 

Vjix)= [ e'<''^>ZF^idX), JG {1,2}, (3.1) 

where EZf^{Ai) = 0, EZF,iAi)'Z^^{A2) = SiFj{Ai n A2), for any j, k G {1,2} and 
Ai, A2 G i3(]R") (Si is the Kronecker symbol). 



7 



By the above spectral representation for the initial data fl2.2p . we can describe the 
vector-solution {{u{t,x),v(t,x)),t > 0, x G M"} by stochastic integration : 

Proposition 1 Lei w(t, Wo(-)) := {u{t,x;uQ{-),v(t,X]Vo{-)), Wo(-) = {uo{-) , Vo{-)) , be 
the vector-solution of the initial value problem h2.1\) h2.2^) . of which satisfies Condition 
A and B, then 

w(t,x; wo(-)) = g(t;rfi,rf,) / e^<A,-e-*l^l'^(^^l^l^)^ ( f ^j^J^^ \ , (3.2) 
with the inter-relative covariance structure 

Ew{t, x; Wo(-))w^(t', x'; Wo(-)) 

= / e^<^'-'>e-(*+*')l^l"(^+l^l^)^Q(t;ci,,rf,)f''^[f^^ ^ Q{t ; d,,d,f , (3.3) 

JR" \ F2{dX) J 

where Q{t; di, 1^2) is defined in l{2.6\) . 

3.2 Subordinated Gaussian initial data 

In this subsection, we assume further that the initials are subordinated fields, as follows: 
Condition C. We consider the random initial data (12. 2p wo(x) := {uo{x),vo{x)) = 
(?7i(x), ?72(x)), X G M", satisfies Condition B and each component has the following form 

r],{x) := h,{Q{x)), X G M", j G {1, 2}. (3.4) 

The Ci(^) C2{x) are independent, mean-square continuous, homogeneous and isotropic 
Gaussian random fields, each is of mean zero and of variance 1, and for each the spectral 
measure Fj{dX) has the (spectral) density /j(A), A G M", and fj{\) is decreasing for 
|A| > Ao for some Aq > and continuous for all A 7^ 0, j G {1,2}, respectively. Moreover, 
we assume that hj{-), j G {1,2}, are real non-random Borel functions satisfy 

Eh%Cjm<oo, JG{1,2}. (3.5) 
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Under Condition C, we have the spectral representations for the sample paths of 
Cj{x), j G {1,2}, as below: 

Qix) = [ e'<^'^>Jf~iX)W,idX), X G M", j G {1,2}, (3.6) 

where Wj{A) is a Gaussian noise measure, Wj{A), A G i3(M"), are centered Gaussian 



with EWi{dX)Wj{diJ,) = 616{X - /i)rfA rf/i. 

Due to fl3.5p in Condition C, we can consider the following orthogonal expansions of 
hj{u) in the Hilbert space L'^{R,p{u)du) with p{u) = 

hAu) = Ci^^ + f:ci^^^, JG{1,2}, (3.7) 

(T=l * 



where 



(3.^ 



C(J)= [ h,{u)^^p{u)du, JG {1,2}, 

and {H„{u), a = 0, 1, 2, ...} are the Hermite polynomials, that is, 

H^{u) = {-ire'^-^e^, for a G {0, 1, 2, ...}. 

It is known that the following two important properties hold ( see, for example. Major 
[25] . Corollary 5.5 and p. 30 ) : 

E[H,,iQiyi))H.ACf{y2))] = Sj's:ia,\Ri;iy,-y2), yi,y2 E M" (3.9) 

and 

H,{Q{x)) = I e^<^''^^-^'^> n Jf~iK)W,{dX.). (3.10) 

In the integral representation (13.101) . the integration J means that it excludes the diag- 
onal hyperplanes Zi = TZj, i,j = I, p,i ^ j. 
The Hermite rank of the functions hj{-) is defined by 

m,- := inf{a > 1 : C^^V 0}, jG {1,2}. 

Specializing Proposition 1 in Subsection 3.1 to the present subordinated Gaussian ini- 
tials, we have 
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Proposition 2 Lei w(t, x; Wo(-)) := {u{t, x;uo{-)),v{t, x;vq{-))), t> x E M"} be the 

vector- solution 5\) of the initial value problem of which satisfies Condition 

A and C, then all the statements in Proposition 1 remain valid, with \3. ^)) is expressed 
as 



w(t,a;; wo(-)) =Q(t;di,d2) 



(1) 



(2) 



peN 



^j<x-,Ai+-+Ap>-pt|Ai+-+Ap|°(l+|Ai+-+Ap|2)T 



Zp^dX) 

zfAdX) 



where 



Zf^{d\) 



^^([^/h{K)Wi{dK)^ 

(7=1 

^{{^h{K)W2{d\„) 



(7=1 



and the coefficient C^p \ j G {1,2} is defined in l{3.8\) . 



We also impose the following assumption which is related to the long-range depen- 
dence of the underlying Gaussian fields Cj(^)) J ^ {1)2}; we refer to [Ulin] for the notion 
and the literatures of long-range dependence. In the following and henceforth, the no- 
tation /(■) ~ g{-) means that the ratio f{ )/g{ ) tends to 1, as the indicated variable "•" 
tends to oo or tends to 0, according to the context. 

Condition D. The Gaussian random fields Cj{^)i j ^ {1; 2}, in Condition C, have their 
covariance functions to be regular varying at infinity in the sense that: 

L{\x\) 



R.{x) 



n 



\x\^ 



as \x\ — )■ oo, < < — , j G {1,2}, 



(3.11) 



where L : (0, oo) — (0, oo) is a slowly varying function at infinity and is bounded on each 
finite interval; recall that L is said to be slowly varying at infinity if limj^_j.oo [L{cy) / L{y)] = 
1 uniformly for any c G (a, 6), < a < 6 < oo. Here, the relation " ~ " in (13. lip means 
that lim \x\'^^ R(.{x)/ L{\x\) = 1; the similar notation will be used in this section and 

\x\—^oo 

also in Section 7. 

Under Condition D, by a Tauberian theorem (see, for example, the book of Leonenko 
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[TBI p. 66]), the spectral density functions of the random fields Cj{^)y 3 ^ {1;2}, are 
regular varying near the origin as follows: 

/,(A)~K(n,«:,)|Ar^-"L(|^|), as A ^ 0, jG{1,2}, (3.12) 
where the Tauberian constant Kin, = — — v^—r—. 

We note that, for each natural number p > 2, the power of the covariance function 
{Rq (x))'' itself is still the covariance function of some random field, for which there exists 
the corresponding spectral density function {fj)*P{X). Indeed, the function (/j)*^(A), A G 
R*^, is the p— th convolution of fj{X) defined as 

{fjT'W = / /.(A - Ai)/,(Ai - A2) ■ - \,-,)a\,-i)lld\i. (3.13) 

for p > 2. Since L''(|x|) is still a slowly varying function for any p, when the p satisfies 
< pKj < n, we can apply the Tauberian theorem again to get 

ifjT^W ~ K{n,pKj)\X\P''^~"LP{\h), as A ^ 0, < pKj < n, (3.14) 

A 

for j G {1,2}. While, if puj > n then the covariance function {R^.{x))'' belongs to the 
class L^(M"); thus the corresponding spectral density function is everywhere continuous 
and satisfies 

(27r)"(/,)*^(0) = [ {Rc,nx)dx< f \R<:,{x)\Pdx < [ \R^^{x)f dx < 00, (3.15) 

Jr" Jr" JR" 

where p* := inf{p G N| pnj > n}; we note that |_R^^. (■)! < 1. 

The displays fl3.13p . (13.141) and (13.151) will be used in the proofs in Section 7. 

4 Micro- scalings for the solution vector-field 

In this section, we present the main result of this paper, which concerns with the micro- 
scaling of the homogenization of the spatial-temporal random field associated with (2.1), 
with the initial data (2.2) subject to the conditions in Section 3. We show that both the 
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Riesz parameter a and the Bessel parameter 7 plays their roles in the scaling procedure. 
The results in this section are more subtle than the macro-scaling discussed in the next 
section; see the remark below Theorem 1 for the interpretation. 

Firstly, we prove the following micro-scaling of homogenization for a single fractional 
kinetic equation, subject to the random initial data. 

= -/i(J- A)i(-A)t,(t,x), .(0,x) = /i(C(x)). (4.1) 

To our knowledge, the homogenization present in the below is a completely new type. 
In (4.2) below, the notation imposed on wants to mean that the variable of C is under 
the indicated dilation factor. 



Theorem 1 Let s := s{t,x; so{-)), t > 0, x E M", be a solution of ( [^. j[ ), which satisfies 
the above Condition B, C and D with k, G (0, ^), where m denotes the Hermite rank of 
the non-random function h{-) on M, which has the Hermite coefficients Ci{h), i = 0,1, . . . 
(i.e., hi{x) = h{x), Ci{^) = C{^)y fiW = fW '^fid ni = n, etc. in Section 3). Then, 
for any fixed parameter x > 0, 

(1) The behaviour of the covariance function of the rescaled random field s^{t,x), t > 
0, X e M", 



s'{t,x) := [e""'''L"'{e~'')]~^\s{et,e—x;h{C{e~— -Co{h) \, (4.2) 
is given by: 

e^<^-^'^> '' dr. (4.3) 



CK + 7 



(2) When e — )■ 0, the rescaled random field s^{t,x), t > 0, x E M", converges to the 
limiting spatial-temporal random field Sm{t,x), t > 0, x G M", in the finite dimensional 
distribution sense, and Sm{t,x) is represented by the Multiple- Wiener integrals 

sUt,x):=^K{n,K)f U^idz,), (4.4) 
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where f ■ ■ ■ denotes a m-fold Wiener integral with respect to the complex Gaussian white 



noise W{-) on . 

Remark. To compare with Proposition 4 in the next Section 5, Theorem 1 has 
the features that it involves both the Riesz and the Bessel parameters, and that it also 
needs to rescale the initial condition. The intuitive meaning behind the latter situation 
is that, while the micro-scaling enforces to "freeze down" both the time t and the space 
X, besides the overall renormalization we also have to "heat up" the initial data, in order 
to get a non-degenerate (though singular) limiting field. 

Now, the micro-scaling of the system is 

Theorem 2 Lei w(t, a;; Wo(-)) := {u{t,x\UQ{-)),v{t,X]VQ{-))), t > 0, x G W, he the 
solution- vector of the initial value problem and Ii2.^) . satisfying the Condition A, 

B, C and D. In the following, x ^■^ o positive parameter, Q{t] di, ^2) is a matrix defined 
in Ii2. 6\) . and the two Gaussian noise fields Wj, j G {1,2} are totally independent. 
Additionally, mi, m2, Hi and K2 denote the parameters contained in Condition G and 
D for uo and vq. 

(1) If m2K,2 > miHi, then the finite- dimensional distributions of the rescaled random 
field 



[£™^i:'"i(£-'^)]-Uw(£t,£^x;wo(£-^-'^-))-Q(^^;^i,^2) J2) |, t>o, xeM", 



converge weakly, as e — >■ 0, to the finite- dimensional distributions of the random field 





where 




with Wi{-) is a complex Gaussian white noise on M" ( i.e., ^.4\) with m, k and W 
replaced by mi, Hi and Wi, respectively ). 
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(2) If rriiKi > m2K2, then the finite- dimensional distributions of the rescaled random 
field 



[e"'''''''L"''{e-^)]-^[w{et,e^x;wo{e-^-''-))-Q{et;d,,d2) ^^2) j }, t>0, xe 
converge weakly, as e ^ 0, to the finite- dimensional distributions of the random field 

Y**{t,x) i ■ [ X^l{t,x) 



t > 0, X e 



where 



and W2{-) is a complex Gaussian white noise on M" ( i.e., Ili4.4\ ) with m, k and W 
replaced by m2, ^2 and W2, respectively ). 

(3) If mi = m2 := m, = K2 '■= k, then the finite- dimensional distributions of the 
rescaled random field 

[e'^''''L"'{e-^T-^{w{et,e^x-wo{e-^-''-)) - Q{et-,di,d2) j }, t > 0, x G R", 

converge weakly, as e ^ 0, to the finite- dimensional distributions of the random field 

Yr*{t,x) Xi^\t,x 
^ Y2***{t,x) J '~ \ Xi\t,x 

where Xm and Xm , are defined in ( [^.5| ) and (^T^ with rrii = m2 = m and Ki = K2 = k. 

To understand the stochastic structure of the hmiting fields, we state, for instance, 
the following covariance result of {Yi**{t,x),Y2**{t,x)). 

Proposition 3 For each fixed t>0, the limiting vector field ^ Y{**{t,x) Y2**{t,x) j 
in the case (3) of Theorem\^is spatial-homogeneous and its covariance matrix has the 
following spectral representation 

^ { ) ( yr{^,^') Y2-it\x') ) = / e^<— '^>5(A;«,7)rfA, 



, t > 0, X G M", (4.7) 



a+7 



where S{X;a,-f) := K{n,mK Y ^ f '* {^(J^'^^f 
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Remark. In view of the singularity of the spectral matrix near the origin, we may 
conclude that, for limiting vector field in the case (3), the long-range dependence (LRD) 
not only exists for each component field but also exists between the two component 
fields; this is a rather new phenomena for LRD, to our knowledge. Similar situation 
happens for other cases, too. 



5 Macro- scalings for the solution vector-field 

In this section, we present the macro-scaling limits for the solution of the fractional 
kinetic systems (12. ip and (12. 2p . in which only the Riesz parameter a plays its role in the 
scaling. 

We again begin with the following single-equation case, which is adapted from [3l The- 
orems 2.2 and 2.3]. 

Proposition 4 Let s := So(-)), t > 0, x G , satisfies ( [^. which satisfies the 

above Condition B, C and D with k, G {0,n/m), where m denotes the Hermite rank of the 
non-random function h{-) on M, which has the Hermite coefficients Ci{h), i = 0,1, . . . 
(i.e., hi{-) = h{-), (i{x) = ({x), fj{X) = /(A) and ki = k, etc.). Then, 

(1) The behaviour of the covariance function of the rescaled random field s^{t,x), t > 
0, X e W, 

s%t,x) :=[e^L™(£-i)]-i[,(i, 4; /,(((.))) -Co(/i)], t>0, xeR^^ 
is given by: 

/ i<x-x' ,\>-iJ,{t+t)\\\°' 
(5-1) 

(2) Moreover, the finite dimensional distributions of the rescaled field s^{t,x) converge 
weakly, as e ^ 0, to the finite- dimensional distributions of the random field 

sUt,x):=^K{n,K)f U^idz,), (5.2) 
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X G M", t > 0, where f ■ ■ ■ denotes a m-fold Wiener integral with respect to the complex 
Gaussian white noise W{-) on M". 

Remark. The above (15.1 p is expressed on the "Fourier-domain " , which is more suitable 
for we will need; while that (2.40) in [3] is in term of the variable domain. 
Then, the macro-scaling of the system is 

Theorems Let w{t,x;wo{-)) := {u{t, x;uo{-)),v{t, x;vo{-))), t > 0, x E W\ be the 
solution-vector of the initial value problem Ii2.1\) and Ii2.2\) . satisfying the Condition A, 
B, C and D. In the following, Q{t; di, 62) is the matrix defined in 1(2. 6]) . pij is the entry 
in Ili2.4\ ), and the two Gaussian noise fields Wj, j G {1,2} are totally independent. 
Additionally, mi, m2, ki and K2 denote the parameters contained in Condition G and 
D for Uq and Vq. 

(1) If ^2^2 > rniHi anddi > d2, then the finite- dimensional distributions of the rescaled 
random field 

[£^L™H^'«r^e-'^^Uw(-,4;wo(-))-g(-;cii,t/2) ( ^ ) j, t > 0, x G M^ 

converge weakly, as e ^ Q, to the finite- dimensional distributions of the random field 

PnP22X^l{t,x) 

^ P2lP22X^nl{t,x) 

where 

r'(l) /•' pi<X,Zl + ... + Zm,>-^J.t\zl + ... + Zm,\" ^ 

'-'»ni f _ \^^^ I ^ TT 

(kl|---kmj)^ 1=1 

with Wi{-) is a complex Gaussian white noise on M" ( i.e., Ii5.2\) with m, k and W 
replaced by mi, ki and Wi, respectively ). 

(2) If mi^i > m2hi2 anddi > d2, then the finite- dimensional distributions of the rescaled 
random field 



(^(1) r' i<X,Zl + ... + Zm-i>-lJLt\zi + ... + Zmi\'^ "!L 

X^^%x):=^K{n,Ki)^ ^ — --—^^ \{Wi{dz,l (5.3) 



[e'^L^-{e-i)]-h-'^^w{i,^■,wo{■)) - Q{-;d„d2) [ ^2) | \, t > 0, xe 

^ e°' ^ \ Gq 



(1) 
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converge weakly, as e — )■ 0, to the finite- dimensional distributions of the random field 

-puPi2Xinl{t,x) 



X^l{t,x):=-^Kin,K2)^ — —^^ l[^2{dzi), (5.4) 



„2 ^' gi<X,2l + ... + 2„2>-Ati|zi + ... + 2„2|" ^ 

anc? Vr2(") a complex Gaussian white noise on M" ( i.e., 115. 2\) with m, k and W 
replaced by m2, H2 o-nd W2, respectively ). 

(3) If mi = m2 := m, /ti = K2 := k, and di > d2, then the finite- dimensional 
distributions of the rescaled random field 

[£^L-(£-^)]-^e-'^^^{w(^,4;wo(-))-Q(^;cii,rf2) j }, t > 0, a; G M^ 

converge weakly, as e 0, to the finite- dimensional distributions of the random field 

T^^\t,x):=T^^\t,x) + Tg\t,x), (5.5) 

where Tm\t,x) and Tm\t,x), are defined in the case (1) and the case (2) with mi = 
m2 = m and ni = k,2 = i^- 

Remark: In the above, we assume that di > In case di < d2, all the correspond- 
ing assertions hold, by interchanging the roles of mi, ki and m2, ^2, etc. As for di = d2, 
it is reduced to the uncoupled case and the result is induced from Proposition 4 directly. 



6 Time-fractional systems 

We extend the above results to the time- fractional derivative /3 G (0, 1), in the 
system (12. ip . that is. 



dt^ 



= -^(/-A)2(-A)2 ]+Bi «, 7>0, (6.1) 



V 
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We recall that the time-fractional derivative ^ is defined (see, for example, the book 
of Djrbashian [H]) by, for any /3 > 0, 

d^f.,,_\ f-'^'Kt) if /3 = mGN 

y r(m-/3) Jo (t-r)/^+i— II P G - i,mj, 
where f^'^\t) denotes the ordinary derivative of order m of a causal function f{t) (i.e., 
/ is vanishing for t < 0). 

The solution of (6.1) can be obtained, under Condition A, by applying the Laplace 
and the Fourier transforms (see, for example, [231 121]), as follows. 

w(t,x;wo(-)) = / P\ ^ ^ \ (\\'^y^ (6-3) 

with the fractional Green function Gi3{t,x; dj) is defined by the transformation 

E^(-/i|A|"(l + |A|2)it'^ + d/) = f e<^^^>Gp{t,x-d,)dx, JG{1,2}, 

where Ep{-) is the Mittag-Leffler function defined by (see, for example, [3] or [21 Chapter 
1]) 

and we shall use the following basic properties about the Mittag-Leffler functions: they 
are entire functions on the complex plane and their asymptotic behavior, when /3 G 
(0, 2], /3 7^ 1, 2, has the inverse power law as follows: 

\Ep^^{z)\ ~ '^(^)' l-^l ^ oo with |arg(-z)| < 7r(l - ^), V 7 > 0, (6.5) 

where arg: C (— tt, tt) and the notation f{z) ~ 0{g{z)) means that f{z)/g{z) remains 
bounded as z approaches the indicated limit point; see, for example, the classic book by 
Erdelyi et.al. [II] (pp. 206-212, in particular p. 206 (7) and p. 210 (21)). 

We note that, from [211 (45)], for /3 G (0, 1), there is an another representation for 
the fractional Green function 

POO 

Gp{t, X] dj) = I fp{t-^s)G{s, X] dj)ds, t>0, xeW, (6.6) 
Jo 



where 

1 n j 

G{s,x]d.) = ( )2e"^e'^^^ s > 0, x eW^, (6.7) 

while f{p),p > 0, is a probabihty density which can be represented by the H-function 
(see, for example, [311 Section 3] and [271 P- 284]) and its Laplace transform is given by 



/■oo 

/ e-'^'Ms)ds = Efsi-q), q > 0. (6.^ 
Jo 



Hence, 



f Gis{t,x;dj)dx^^ j t ^ j ffsit ^ s)G{s,x]dj)ds dx 



oo 



=t ^ I fisit ^s) I G(s, (ij)(ia; (is, (by Tonelli theorem) 

E3 



1 ("Xl \f? ol 

^t-^ / Mt-^s)e'^^ ds ^ Ep{d,tP), 
Jo 



(6.9) 



C(t;i?):=i?w(t,x;wo(-)) = P| ' ^ , , \P'\ ^% |, (6-10) 



where the convergence of the integral in (16.91) is guaranteed by the asymptotic behavior 
of the H-function (see, for example, [271 (3-7)]). 

From the above discussion we know G/3(t, x; dj) G L-'^(]R", dx) and Gjslt, x; dj)dx = 
Eji{djtl^) for any t > and j G {1,2}. Therefore, if the initial data Uq{x) and Vq{x) 
satisfy the special form (13.41) . then by the representation (16.31) we have 

Ep{d,t^) \ / Cj^) 
Ep{d2t^) ) \ 

where j G {1,2} are the Hermite coefficients defined in fl3.8l) . 

The following Theorems 4 and 5 are time-fractional versions those in Theorems 3 and 
2, respectively. However, it is now needed to add the feature of the sub-diffusive property, 
which is the reflection of time-fractional /3 < 1 (see Section 1), into consideration of the 
macro-scaling of homogenization. We need to take an additional scaling on the matrix 
B in the system (16. ip in order to compromise the effect of this sub-diffusivity upon 
the interaction between u and v. To emphasize this situation, we denote the vector 
solution w(t, x\ Wo(-)) by w(t, x; Wo(-), B) in the following formulation of macro-scaling 
of homogenization of a spatial-temporal fractional kinetic system. 
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Theorem 4 Let {w(t, x; Wo(-), 5), t > 0, x E M"} he the solution-vector of the initial 
value problem \6.1\) and satisfying Condition A, B, C and D. Moreover, the LRD 

parameter kj and the Hermite rank rrij satisfy rrijUj < mm{2a,n} and the Gaussian 
noise fields VTq) are totally independent for j G {1,2}. 

(1) If rriiKi < 7712^2, then the finite- dimensional distributions of the rescaled random 
field 

t > 0, X E M", converge weakly, as e — 0, to the finite- dimensional distributions of the 
random field 



x) = I ^11^22^(1) (t, x; di) -pi2P2iT(^)(t,x;d2) 
\ P2iP22T^^^ {t, x; di) - p2iP22T(^) {t, x; 62) 



t > 0, X e 



where for j G {1, 2} 



T^^\t,x;dj) (6.11) 

(2) // 7722^2 < rriiKi, then the finite- dimensional distributions of the rescaled random 
field 

T^^\t, x){e'^ L"^' (e-^))'^ \^w{e-\ x; wo{-), B) -C{e-H; e^B)^ 

t > 0, X G M"", converge weakly, as e ^ 0, to the finite- dimensional distributions of the 
random field 

\ -pi2P2iT^ ' [t, x;di)+ puP22T'^ ' [t, x;d2) J 
where for j G {1, 2} 

T^^\t,x;dj) (6.12) 
Vm2^- J flAil---IA^.n^ 



m2\) " 1=1 
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(3) If nil = 7712 = m and Ki = K2 = n, then the finite- dimensional distributions of the 
resettled random field 

Tf)(t,x) := {e'^V^{e-^))-'^\w{e-H,e-^x-WQ{-),e^B)-C{e-H-e^B)y 

t > 0, X G M"", converge weakly, as e — > 0, to the finite- dimensional distributions of the 
random field 

T^'^\t,x) =f^^\t,x)+f^^\t,x), t>0, xeM", 

where the random field T^^\t,x) is the same as the limiting random field T^^\t,x) in 
( I6'.ii)) and 116.1^) by replacing rrij — )• m and Kj k for j G {1,2}. 

Remark. The restriction mjKj < min{2a, n} in the above Theorem 4, together with the 
power law decay of Mittag-Leffler functions, guarantee that the random fields T^^\t, x), j E 
{1,2,3}, are indeed defined as L'^{Q, J^,V) stochastic integrals. 

As for the micro-scaling, the sub-diffusivity has no influence, and the same micro- 
scaling procedure as Theorem 2 applies. 

Theorem 5 Let {w(t, x; wo(-), 5), t > 0, x G M"} be the solution-vector of the initial 
value problem 116. 1\} and 112. satisfying Condition A, B, C and D. Moreover, the LRD 
parameter Kj and the Hermite ranknij satisfy rrijUj < min{2(a-|-7), n} and the Gaussian 
noise fields W(j^ are totally independent for j G {1,2}. 

(1) If rriiKi < m2K2, then the finite- dimensional distributions of the rescaled random 
field 

^(1) ._ wo(5"^"^-)) " ( ^(2) 1 }> ^ > 0, a; G M", 

converge weakly, as e 0, to the finite- dimensional distributions of the random field 
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where 



J (|Ai|---|A„J)^ l\ 

Jgnxmi VI J- I I "ill/ ( — 1 

(2) // 7712^2 < rriiKi, then the finite- dimensional distributions of the rescaled random 
field 

( r^^^ \ 

:= [£™2''2^L'"^(e">^)]-5|w(et,e^a;; wo(£"^"^-)) " ( ^(2) j }> ^ > 0, x G M", 
converge weakly, as e ^ 0, to the finite- dimensional distributions of the random field 

M(^)(t,x)=| l,t>0, XGM^ 



where 



(3) If mi = m2 = m and Ki = K2 = 1^, then the finite- dimensional distributions of the 
rescaled random field 

( r^^^ \ 

M(3) := [£-->^L"^(£-^)]-i|w(£t,£^x; Wo(£-^-'^-)) " ( c'f ) }' ^ ^ °' ^ 
converge weakly, as e — 0, to the finite- dimensional distributions of the random field 

where the random field M(-'^(t, x) is the same as the limiting random field M*^''^(t, x) of 
Case (1) and (2) by replacing nij — )■ m and — a for j G {1,2}. 

The concluding remark: The time-fractional index /3 < 1 indicates the sub-diffusivity, 
and it changes to be the super-diffusivity if we consider /3 > 1 (see Section 1). In [22], 
the time-fractional reaction-wave type system with random initial data are studied, in 
which the first-order time-derivatives of the initial data play the crucial role. To consider 
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spatial-temporal fractional kinetic systems which is super-diffusive in time and Riesz- 
Bessel in space will be a task of tremendous analysis. Finally, we mention that, for 
the classical, i.e. non-fractional, heat-type system with random initial condition, the 
solution vector-field and the scaling limit are expressed in terms of heat kernels; this 
more explicit and simpler case is treated in pT| . 



7 Proofs 

In the following proofs, ^ denotes the convergence of random variables (or random 
families) in distributional sense, and = denotes the equality of random variables (or 
random families) in distributional sense. Moreover, we also denote f{t,x]e) x g[t,x]e) 
if there exists a constant c := c{t,x) > such that cg{t,x;e) < f{t,x;e) < c~^g(t,x]e) 
when e — )■ 0. 

Proof of Proposition 1. 

For (13. 2p . we use the solution form (12. 5 p and Karhunen's representation (13. ip to get 
wit,x;w,{.))=Q{t,d,,d,) [ G{t,y;a,^) ( '7<TZll^'^!!i^, I ^^'^^ 
= Q{t, du d,) [ [ e^<''^-y>G{t, y- a, ^)dy ( [ ) 

JR" \ ZF^dX) J 
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For fl3.3p . it is a consequence by using the independence assumption between the initial 
data, and we proceed it as follows, 



Ew{t, x; Wo(-))w(t', x'; Wo(-)) 

--Q{t; di, d2) [ e*«A,->-<A',x'»e-/.t|Ar{i+|AP)5V|A>(i+|A'P)5 ^ 
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Before going to prove our main results we recall the following two arguments, which 
are powerful to help us to reduce and simplify our problems. 

(Slutsky argument) (see, for example, the book of Leonenko [161 P- 6.]) Let {^e} and 
{?7e} be families of random variables such that {^s} =^ ^ and {?7e} c, where c E R. 
Then {i) ie + Ve ^ i + c, {ii) i^Ve c^, and {in) ie/Ve i/c 

(Cramer- Wold argument) (see, for example, again [111 p. 6.]) A family of k- 
dimensional r.v's := (^e(a;i), ^^(xfc))^ converge in distribution to a /c- dimensional 
r.v. ^ := (^(xi), ...,^{xk))'^ if and only if, for every c := (ci, ...CkY ^ "^^^ 

k k 
< C, >= ^ Ci^e{Xi) => ^ Q^(Xi) =< C, ^ >, aS £ 0. 
i=l i=l 

The following lemma, although it is a simple extension of Cramer- Wold argument, 
is of interest in itself, and will play an important role in the proof of our Theorem 2 and 
3. 

Lemma 1. Let := {[xP , xPf{t,x), x G M", t > 0} be a M^.yalued random field 
which is generated by X^(t,x) = Qe(t)[t/e, V^]"^(t,x), where U^it^x) and V^(t, x) are 
independent random fields on R"" x M'^ and Qe{t) is a non-random 2x2 matrix. If there 
exist two random fields Uq and Vq such that Ue{t, x) =^ Uo(t, x) and Vs{t, x) =^ Vo(t, x), 
respectively, and Qeit) converges to Q(^t) in the usual sense when e — )■ 0, then the finite 
dimensional distributions of Xi;{t,x), t > 0, x E M", converge to the finite dimensional 
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distributions of X := {Q{t)[Uo, Vof{t,x), t > 0, x e W}. 
Proof of Lemma 1 

By Cramer- Wold argument with k = 2 there, it suffices to prove: For any given Ci, 
C2 G M and x, t fixed, we have 

[ci, C2][g.(i)][f/.(x,t), K(x,t)]^^[ci, C2][Q(t)][Uo{x,t), V,{x,t)Y, 

which is equivalent to 

(ClQll(t) + C2Q2l{t)) Uo + (CiQuit) + C2Q22{t)) Vo, 

where the i,j indicate the entry of the matrix. While the above display can 

be checked by using the characteristic functions, since f4 and are assumed to be 
independent (whence so are [/q, Vq). □ 
Proof of Theorem 1. 

(1) Firstly, for simplification, we set Ni{e) := e^"^'^L"^{e~^). By the Hermite expansion, 
we can rewrite s{et, x; h{({e~"+^~^-))) — Co(^) as 

s{et, e^'x; /i(C(£"^"^-))) " Co{h) = V s{et, e^x; ^^i7JC(e"^"^-))), (7-2) 

VP- 

p=m * ' 

where the summation is in L'^{VL) sense. Hence, in accordion to the definition fl4.2p about 
the random field x), it can be rewritten as 

oo 

s^(t,x) = 5^/;(t,a;), (7.3) 

p=m 

With j;(t,x) := {N^{e)y-2s{et,e^^x- {({8-^ -''■))). 
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By fl2.5p with rfi = ^2 = 0, we have 
'^'^4^-\m{e)f^i;{t,x) (7.4) 

G{et, y; a, -f)Hp{C{e~^~^{e^ x - y)))dy 



(T=l 



1 



For the bracket above, by substituting t ^ et and A — )■ e "+i ^{Xi + ■ ■ ■ + Xp) into (12.31 
we have 



G{et,y;a,'y)e 
so (I7.4p is equal to 



^i<e-Xx,Ai+-+Ap>g-Mete ^ "^|Ai+-+Ap|"(l+e ^ '''|Ai+-+Ap|2)7 "Q ^ J(;^^)p^(^_X^j 



(T=l 



' J eK.x+---+A;>e-/..te-"+-iA;+...+A;r(i+.-4.|A;+...+A;p)5 j^^^ 

(7.5) 

where we have used the self-similar property for Gaussian random measure on in the 
last equality. Therefore, by the orthogonal property for the Gaussian white noise, we 
can get 



iCpih)r'N,ie)Covimx)F(t\x')) 



, a 2 7 P 

rXpn I i<x-x ,X^+-+Xp> -^le{t+t)e "+T |Ai+---+Ap|"(l+£ °+t |Ai+---+Ap|2)7 



\[f{e''X,)dX, 



(7=1 



Xpn I i<x-x' ,Ti> ~iie{t+t)e ^|ri|°(l+e ^|ri|2)i 



5X(P-1)'^ 



2 T 
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where f *''{■) is defined in fl3.13p . 

(i) For p G N with mn < pn < n and any 5 > 0, by fl7.6l) and fl3.14p . 

Cov(j;(t,x)j;(t',x')) = {C,{h)f{A,{e) + A^{e)), (7.7) 

with 

< {N^{e))-h''^ sup{r''(A)| |A| >5] f e-^(*+*')l"l"^'dr 

J\£^t\>5 

— )■ 0, as £ — )■ 0, {niK < pK <n) 
and, by choosing 5 small enough and ( I3.14p . 

J\e^T\<5 

= (iVi(£))-V^ /" gKx-x',.>g-MW ^ o(l))/^(n, p/c) ^ '^^jj",'^ dr 

|eXr|<<5 

/ i<X-x' ,T>-ll{t+t)\T\°'+'^ 
as £ 0, 

where the asymptotic equivalence is guaranteed by the uniform convergence theorem for 
the slowly varying function (see, for example, Leonenko [121 Section 1.4]). We remark 
that the /, defined as a spectral density function, is bounded outside of zero(while the 
singularity at zero is from the LRD assumption, as employed in subsection 3.2); therefore 
its p-th convolution remains to have, at most, the singularity only at zero. 
The conclusion of (i): Apart from the term Cov(/^(t, x)/^(t', a;')), 

lim Cov{i;{t,x)i;{t,x)) = (7.8) 

p:m<p<n/K 
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where we have used the fact that {/ G N|m/t < Ik < n} is a finite set and on this set 
hm(Ni(e))-^e^p''LP(e-^) = except for p = m. 

oo 

(ii) For p G N with > n , hj fIXTSD and ^ {Cp{h)f <\\ h g< oo, 



hm J2 CoYii;{t,x)i;{t,x)) (7.9) 

p:pK>n 

=hm e"'^(iVi(e))^i ^ {Cp{h)f j e^<---''->"A'(*+*')l-l'^+7*^(0)dr 

p:pn>n ^„ 

<hm e"'^(iVi(e))-^M V (C„(/i))2 = 0, 



p:pK>n 



since by fl3.15p /*^(0) is bounded by /*''(0) with p = mi{l e N\ Ik > n} so we set 

M:= /*?(0)/ig„e-'^(*+*')H"^"rfr. 

Finally, from the expansion (17.31) for the random field s^(t, x) and combining the obser- 
vations (17.81) and (17.91) we know that only the component I^(t,x) in (17.31) do contribute 
to the covariance function of the random field s'^{t,x), that is. 



Ihn Cov{s%t,x)s%t,x')) = {CUh)?K{n,mK) I e'<^-^' '->^-—^dT. □ 



-lM{t+t')\T\°' + 1 



(2) From the above discussion, we may apply Chebyshev inequality to obtain that: 

oo 
p=m+l 

Therefore, in view of Slutsky argument, we suffice to focus our attention on the term 
/^(t, x). In the following we will prove /^(t, x) converges in distribution sense to Smit, x), 
which is defined in (14.41) . for each fixed (t, x) G M+ x M". By the definition of Ni[e) and 
replacing the letter p by m in (17. 5p . we can rewrite (I7.5p as follows 

Ut.x) = / e^<^'<+-+^''">M,(A) ni^«), (7.10) 



with 

2 

M,(A) := e^^^L-f (e-^)e-^*(^^+l^^+-+^-l')^l^^+-+^-l'^ J] ^/f^) 

CT=1 
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which, when e — > 0, satisfies 

hmM.(A) ^ {K{n,K))f ' (7.11) 

(|Ai|---|A^|) — 

Now, applying the isometric property of the muhiple Wiener integrals to the difference 
of frrrrm and (gai), we have 

E|/^(t,x)-.,^(t,x)|2=lim(Cp(/i))2 / |M,(A)-(K(n,«:))f-— — — ^^-pJJrfA. 

X m 

— y 0, as e — !■ 0, by (17. lip and the assumption /(A) is decreasing at infinity in Condition 
C and 

„-2Ajt|Ai+-A„|"+T ™ /•„-2Mt|A|"+'>' 

(T=l 



r ^-2iit\Xi+-Xm\°'+-' ™ /•g-2Mt|A|"+'^ 

/ n — IT dXa = r(n, m, k) / ,, , ,, dX < oo, for mK < n, 

y (|Ai|---|A„|)-« -Li ' V (|A|)-— 



where the constant r(n,m, k) is generated by the Riesz potential. Finally, the assertion 
(2) of Theorem 1 is followed from Slutsky and Cramer- Wold arguments. □ 
Proof of Theorem [2] for the case (1): m2a2 > niiai. 
From the solution form (12. 5p . we have 



Vit^x) i I Cf 



(7.12) 



where Q(t;di,d2), U(t,x) and V(t,x) are defined in (12.61) and (12. 7p . 

By ciai, 



(1) 



(2) 




if U{et,e°'+-ix;uQ{e "+t ^■)) — Cq^^ 
V{et, e~x; uo{€~~~^-)) — Cq 



Q(.M„d,)[.™L-(^-rM :r: 'T-^-.ii ;(2) 1 (7-13) 
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Firstly, by Theorem [T] (2), we have 

U,{t,x)^X(^l{t,x), (7.15) 

where Xml is defined in fl4.5p . 

Secondly, by Theorem [1] (1), we can obtain 

Veit,x)^0 (7.16) 

since we can apply Chebyshev inequality to observe that for any c > 0, as e — )■ 0, 

P(|K(t,x)| > c) < c-Var(V,(t,x)) x c-2[5-"*i"«L-'"i (e"^)] ■ [^-'"^"^^^-'"^(e-^)] ^ 0. 

Meanwhile, since 

lim Q,{t) := lim Q{et; d^,d2) = p| I M P^' = h^2. (7.17) 

Therefore, we may apply Lemma 1 to those Ui;{t,x), V^{t,x) and Qe{t) on the above to 
obtain that 



, -J- , 1 V ^^ -Wl", "l,"2j (2) 



Proof of Theorem [2] for the case (2): rriiai > m2a2- 

The proof is proceeded as the case (1), yet under the new assumption and the different 
renormalization [e'^^'^^^ L"^^ {e~'^)]~^ . Now fl7.15p becomes as 

Us{t,x) := [£™2'=2^L'"2(e-x)]-i|t/(et^£^a;;^Q(^-^-x.)) 0, 
and fl7.16p becomes as 

V,{t,x) := [£™>^L'"2(£-^)]-i|l^(£t,e^x;i;o(£"^"''-)) -C'^} 
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- — ' /oN 

where Xmi is defined in (14. 6p . 

While f l7.17p is kept unchanged. Therefore, we again apply Lemma 1 to get 



(1) 



v{et,e°'+-t x;vQ{e "+t ^■)) I \ C, 



(2) 








/2X2 I ~(2),, , ], t>0, XeW. □ 



Proof of Theorem [2] for the case (3): uii = m2 = m, ai = a2 = a. 

By Theorem [T] (2), we have 

:= [e"'-^L"^{e~^)r'^{v{et,e^x;vo{e-^-''-))-Ci,'^}^X^\t,x), 

where J G {1,2}, are defined in (14. 5p and (14.61) with mi = m2 = m. 

Because in this case the equality limQi^U) = / is still unchange, in the same way, we 

e-s>0 



obtained 



1 j / u{et, e°'+-' x;uQ{e "+t ^■)) \ j j \ / ^( 



(1) 



(2) 




Proof of Proposition [31. 



Y***{t\x') Y***{t',x') ) 

EX^')(t,x)xi')(t',x'). 

Because the representation for the limiting fields Xm(t,x) and Xm\t,x) is the same as 
the limiting field Sm(t,x), defined in (14. 4p . we can apply the result (14. 3 p to get 

EX(jXt,x)X^i\t\x') = {C^£fK{n,Km) / e^<— ^^>————dT. 



(7.18) 
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Therefore, the covariance structure (I7.18P is equal to 





Proof of Proposition 3] 

(1) Here, for simphfication, we set G{t,x) := G{t,x;a,'y) and N{e) := L'^{e~^) , 
then by the solution form (12 .Tp for the differential equation = — — A)2 (— A)^s 
we have 

iN{e))-Hs{lA;HC{-)))-c,{h)} 

= (iV(5))-^{ / GiK^-y)\Coih) + f2c,ih)^^^^]dy-Coih)] 
Y.'^{N{e))-U G{lA-y)H,{ay))dy:=f2sl{t,x). (7.19) 



fc=m fe=m 

From (13. 9p . the cross terms of the left hand side blow have zero covariance, thus we have 



CoYiJ2'l{t,x),J24it,x')) = J2Coy{sl{t,x),sl{t,x')). (7.20) 

k=m k=m k=m 

For each k G {m, m + 1, . . .}, by (13. 9|) 

Cov(4(t,x),s^(t',x')) 
= (C,(/^))2(iV(e))-i / / G(-, 4 - y)G{-, ^ - y)R\y - y)dy dy 

R" R" 

= {Gu{h)f{N{e))-^j I G(^, ^ - y)G(^, ^-y'){l e''''-''''^r\\)d\]dy dy' 

R" R" R" 

= iC,{h))\N{e))-'l { 1 e'<^'^>G(^, ^ - y)dy 1 e^<-^''^>G(^, ^ - y')rf2/'}r'(A)rfA 
^(C,(/.))2(iV(e))-^ y"e^^Tr''^-/^¥lA|"(i+|AP)^;*^'(A)rfA 

R" 

=(Cfc(/i))2(iV(£))-i / e^<^-^''^>e-^(*+*')l^l"(i+l^^^l')^£^r'=(£^A)rfA, (7.21) 
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by rescaling A into e^X. 

For k : mn < kn < n and any 5 > 0, by fl7.2ip 

Cov(4(t,x),4(t',x')) = {C,{h))'{A^{e) + A2{e)), (7.22) 

with 

< {N{e))-'e^ sup{r'=(A)| |A| > 6} [ ^ e'^^^^+^'^l^l^rfA 

J\ei\\>S 

POO , 

<C{S)e'^L"'{e-^) / ^ e-'^(*+* ^^V^'Mr ^ 0, as e ^ 0, 
and by choosing 5 small enough 

\eaX\<5 



J |e«A|"-'='' 

|ea A|<(5 

/•pj<a;-x',A>-^j(t+t')|A|" 

~ (iV(£))-^e^L'=(e-^)(l + ?)if(n, / — — d\ as e ^ 0, 



|A| 



where the asymptotic equivalence is guaranteed by the uniform convergence theorem 
for the slowly varying function (see, for example, Section 1.4]). So from the above 
discussions in fl7.22p we can conclude that 

^i<x-x\\>~n(t+t )\X\°' 



/ i<x-x ,X>~iJ.(t+t J|A|" 
\X\^. 

(7.23) 



since {/ G N|mfi: < Ik < n} is a finite set and lim(A^(£:)) ^e'''^ L''{e «) = for A; G {/ G 



N|mK < Ik < n} except for the term k = m. 

33 



For k:kK>n,hy fl?:^ and flXT^ 



k:kK>n 



lim J2 CoY{sUt,x),sl{t,x')) = (7.24) 
\ime^(N(e))-\l + o(l)) V (Ck(h)? f e'<--' ^>-^^(t+tM- f*krQ^^x 

£^0 ^ ^ J 

lim5S(iV(e))-i(l + o(l))M ^ {C,{h)r = 0, 



k:kK>n 



since by fl3.15p /*''(0) is bounded by /*''(0) with k = inf{/ G N| > n} so we set 
M := /*^(0) Jjg„ e"^*^*"^* ^l^l^rfA. The proof of Proposition |4](1) is completed by combining 
dlSni), and (17:2411 to obtain 

/ i<x-x' ,\>-fi{t+t')\X\°' 
\X\^^ ° 

(2) This is derived from (1) by the same way as that in the proof of [31 Theorems 2.2 
and 2.3], and thus we omit it. □ 
Proof of Theorem [3] for the case (1): 7722^2 > niiKi and di > d2- 
From fl7.12p . we have 



= e-*iQ(i;rf.,d.)[.-i?^L"«(.-i)l-i ( '^^^^ l^, ) ■ (T.25) 

Firstly, by Proposition H] (2), we have 

U,{t,x) := [e'^L-^{e-^rUu{lA) -Cl,'n ^ XW(t,x), (7.26) 

where Xml is defined in fl5.3p . 
Secondly, by Proposition H] (1), we have 

K(t,x) := [e'^L"^^{e-^r'^{v{-A)-Ci,'n ^0, (7.27) 
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since we can apply Chebyshev inequality to observe that, for any c > 0, as e — )■ 0, 



P{\V,{t,x)\ >c)< c-2Var(K(t,x)) x c-^[e-^ L-'^'ie")] ■ [e^L^'{e")] 0. 
Meanwhile, since di > d2, 



t J ^ t „ / e""!^ 



limQ,(t): = lime-*.Q(-;d,4)=lime-..P(^ e'-'i j 

= p(^ " V-' = f """" """"" 1 . (7.28) 

\ y y P21P22 -P12P21 J 

Therefore, by the independence between Us{t,x) and Vs{t,x) we may apply Lemma 1 to 
obtain 







P21P22 y \ 

PuP22X^l{t,x] 



t > 0, X e R". □ 



Proof of Theorem [3] for the case (2): 1112^2 < miHi and di > d2- 

We use the same scheme as in the proof of the case (1) with the roles of U^it^x) and 
Veit^x) being replaced as follows 



U,{t,x) := [e^L-^(e"i)]-^{f/(-,4)-Cj^n Ao. (7.29) 

and 

:= [e"^L-(e-irU\/(^,4)-Cf I (7.30) 



(1'] 

where Xml is defined in (15. 4p . Additionally, in this case the limit of the matrix Qe{t) 



35 



coincides with f l7.28p so 

P11P22 -puPu \ ( 

P21P22 -P12P2I J \ Xinl{t,X 

-pnPuX^l{t,x 
-PuP2.X^lit,x: 



t>0, xe M". □ 



Proof of Theorem [3] for the case (3): mi = m2 = m, ki = K2 = k, and di > ^2 

By Proposition m (2), we have 

. t X 



U,{t,x) := [e'^L-ie~-^rUu{-A) - Cl,'^] ^ X^J,\t,x), (7.31) 
and 



:= [e'^L-{e-^^)]-'^{v{'-A)-Ci,'^]^X^\t,x), (7.32) 



, t X 

'(1) .^A V(2) 



where Xm and Xm is defined in (15. 3p and (15. 4p with mi = m2 = m. 
Additionally, in this case the matrix Qe{t) is also unchanged so by applying Lemma 1 
to Ue(t,x), Vs(t,x) and Qeit) which are given in (I7.3ip . (17.320 and (I7.28p . respectively, 
we see that the finite dimensional distributions of the rescaled random field 



P11P22 -PllPu \ ( Xlr}{t,x) 

P21P22 -P12P21 I \ X^i\t,x) 



PuP22Xm {t, X) - PllPuXinl {t, x) 
P2lP22^m^ (t, X) - Pl2P2lXi^l {t, x) 



t > 0, X e M". □ 



Proofs of Theorems 4 and 5. 

The proofs can be proceeded parallel to the proofs of Theorem 3 and 2, respectively, 
and thus we leave them to the reader. □ 
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